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Users  of  orbital  elements  are  often  puzzled  by  the  different  ways  in  which 
mean  elements  are  defined.  This  applies  to  all  the  standard  elements,  but  the 
difficulty  is  most  pronounced  in  regard  to  semi-major  axis  and  mean  motion. 

This  paper  relates  mean  elements  to  the  perturbations  of  a  satellite  caused 
by  the  Earth's  zonal  harmonic  J2  ,  and  shows  how  elements  can  be  defined  so  that 
perturbations  in  spherical  coordinates  take  an  especially  simple  form.  The  paper 
concentrates  on  the  first-order  effects  of  Jj  ,  with  formulae  given  to  connect 
the  preferred  elements  to  the  elements  of  Kozai,  but  it  also  considers  the 
extension  to  second  order  in  J2  and  first  order  in  the  general  zonal  harmonic. 
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ABSTRACT 

User*  of  orbital  element a  ore  often  puzzled  by  the  different  vaye  in  which  mean  elements  are 
defined.  This  applies  to  ell  the  standard  elements,  but  the  difficulty  Is  most  pronounced  In 
regard  to  seal -major  axis  and  mean  motion. 

This  paper  relates  mean  elements  to  the  perturbations  of  a  satellite  caused  by  the  Earth's 
zonal  harmonic  J2  •  and  shows  hov  elements  can  be  defined  so  that  perturbations  In  spherical 
coordinates  take  an  especially  simple  form.  The  paper  concentrates  on  the  first-order 
effects  of  J2  ,  with  formulae  given  to  connect  the  preferred  elements  to  the  elements  of 
Kozai,  but  it  also  considers  the  extension  to  second  order  In  J2  and  first  order  In  the 
general  zonal  harmonic. 

1  INTRODUCTION 

It  Is  a  notorious  fact  that  the  orbital  elements  derived  by  the  computer  programs  of 
different  Institutions  are  likely  to  be  defined  In  different  ways,  and  this  Is  a  frequent 
cause  of  confusion  to  users.  The  difficulty  Is  an  obvious  one  when  entirely  different 
element  sets  are  used,  but  It  Is  just  as  serious  when  the  element  sets  are  In  principle  the 
same.  In  the  present  paper  we  assume  the  most  commonly  used  'standard*  set,  In  which 
a,  e,  1,  Q,  w  and  M  are  supplemented  by  n  ,  the  mean  motion.  If  osculating  elements  were 
published,  there  would  be  no  difficulty,  n  not  even  being  necessary  if  the  value  of 
U  (  -  GM)  Is  known;  In  practice,  however,  it  Is  always  'mean'  elements  that  are  published  and 
It  Is  the  precise  meaning  of  these  that  causes  the  problem. 

If  C  Is  any  member  of  the  element  sat,  the  relation  between  osculating  C  and  mean  C  Is, 

C  -  l  +  »  (1) 

where  6c  symbolizes  the  combination  of  all  short-period  perturbations,  whilst  ;  has  long¬ 
term  variation  only.  „lt  Is  immediately  apparent  that  the  spilt  Is  not  uniquely  defined,  _ 
since  the  origin  of  C  Is  essentially  arbitrary  -  thus  any  quasi-constant  component  of  C  , 
in  a  given  split,  can  be  transferred  from  C  to  6c  without  changing  the  nature  of  the 
split,  so  long  as  the  order  of  magnitude  of  this  component  Is  no  greater  than  that  Inherent 
In  6c  . 

We  are  supposing  here  ehat  the  ■  •  rturbations  are  due  to_the  zonal  harmonics,  Ji  ,  and 
will  be  mainly  concerned  with  K  .  .h<n  the  variation  of  C  Is  assumed  to  be  made  up  of 
the  usual  aecular  and  loog-peric^  <  -  Togcthar  with  tha  terms  of  6c  ,  these  can  be 
obcalned  by  Integration  of  Lagrang.  lanetary  equations.  The  variation  of  the  fast 
element,  M  ,  is  obtained  via  introduction  of  tha  subsidiary  element,  0  ,  defined  so  that 

t 

M  •  j  n  dc  +  0  and  M  •  n  +  0  .  (2) 

0 

Direct  integration  of  Lagrange's  equations,  with  t  (or  la  practice  H)  aa  lotegrstlon 
variable.  Involves  the  truncation  of  power  aeries  in  e  ;  if  v  (true  anomaly)  Is  used,  on 
the  other  hand,  a  complete  (untruncated)  solution  can  be  derived  for  the  first-order 
variation  of  C  due  to  Jt  ,  and  likewise  /l /  for  the  second-order  variation  due  to  J2  . 

But  these  two  different  techniques  for  the  Integration  Immediately  lead  to  alternative 
definitions  of  C  ,  on  the  basis  of  s  hypothetical  further  Integration  with  respect  to  M 
and  v  respectively;  If  this  integration  is  taken  over  a  complete  revolution.  It  leads  to 


Che  concept  of  a  *8-m«en^  element  (where  6  la  M  or  v)  such  Chet  the  8-mean  value  of  4;  , 
relative  to  the  8-mean  c  ,  la  aero. 


The  uee  of  H-maan  eleaent a,  aa  In  111  and  / 3/,  la  euperf icially  natural,  but  the  complexity 
and  the  Inaccuracy  aaaoclated  with  the  e-truncation  are  unattractive.  It  would  be  poaalble 
juet  to  adopt  v-maan  eleaenta  lnetead,  but  two  further  conalderatlona  enter  the  picture. 
First,  there  sight  be  an  advantage  If  a  particular  derived  quantity,  viz  r  (the  radial 
dlatance),  had  aa  small  a  blaa  aa  poaalble  relative  to  lta  value,  ?  ,  obtained  from  the  c  ; 
thla  wae  evidently  a  factor  in  the  definition  of  the  widely  uaed  eleaenta  of  Kozal  /A/. 
Second,  a  definition  of  the  c  would  be  eapedally  attractive  if  the  aaaoclated  4c  could 
be  combined  into  a  triple  of  very  compact  expressions  for  perturbation#  in  coordinates.  The 
latter  consideration  motivated  the  study  of  the  present  author  /l/,  based  on  a  particular 
ay stem  of  aphar leal-polar  coordinates  end  on  earlier  work  /5,6/. 

In  the  current  paper  the  uae  of  c  will  normally  refer  to  the  mean  elements  of  / 1/,  so  that 

1  will  be  Identical  with  the  quantity  a'  to  be  Introduced  in  the  next  section.  These  C 
will  be  collated  with  the  mean  elements  of  Kozal,  denoted  by  Cg  for  distinction,  since 
the  wide  ueage  of  the  latter  Includes  their  use  in  the  RAE  orblt-determlnatlon  program 
PROP  /7/. 

2  MEAN  SEMI-MAJOR  AXIS  AND  MEAN  MEAN  MOTION 
2  3 

Since  n  a  ■  u  ,  it  might  seem  that  the  definition  of  S  should  be  automatically  tied  to 
>  ,  being  such  that  .  There  ere  advantages  in  making  II  and  1  Independent, 

2_3 

however,  with  H  a  equal  to  a  modified  w  .  This  was  recognized  by  Kozal  /4/,  whose 
version  of  Kepler* s  third  law  may  be  written  aa  H2!3  ■  u(l  -  K^h),  where 

K  -  |  .  p  -  »q2  -  .<1  -  «2)  ,  h  -  1  -  |  f  .  f  -  .lo2l  ,  (3) 

and  R  ia  the  Earth's  equatorial  radius.  Ironically,  though,  we  shall  find  that  to  first 
order  In  J2  the  preferred  definitions  of  h  and  Z  satisfy  the  unmodified  law;  taken  to 
second  order  in  J2  ,  and  with  the  general  even  harmonic  allowed  for  to  first  order,  the 
modified  law  takes  the  form  (summation  for  even  t  >  4) 

ff2*3  -  pjl  -  *V(8  -  8f  -  5f2)  +  2  VT)  VS)j  >  (4) 


where  the  constants  (Y^ )  and  functions  (A^  and  Q^)  will  be  specified  later. 

In  seeking  a  compact  formula  for  Jr  ,  r  being  one  of  the  spherical  coordinates  referred  to 
in  the  Introduction,  we  have  to  deal  with  the  term  (3r/3a)  4a  ,  where 


3r  r  _ q 

3a  a  l  +  e  cos  v 


(since 


1  +  e  cos  v) 


(5) 


thus  we  would  like  1  +1  cos  V  to  be  a  factor  of  4  a  ,  to  cancel  the  denominator  in  (5). 
Following  the  earlier  work  /5,6/,  we  note  that  conservation  of  energy  (for  the  tonal 
potential  field)  yields  an  exact  constant  a*  (to  whatever  order  the  analysis  is  takan), 
where 

V"  ’  T  +  F'  ’  Vh*r*  U  *  £Ul  *nd  Ut  *  Pl(,In  6)  •  (6> 


U  being  the  disturbing  potential,  8  the  latitude,  a  ia  ft  the  ueual  Legendre  polynomial. 

It  follows  that  a  -  a*  ■  2aa'U/u  and  hence  contains  (R/r)2  as  a  factor  (from  every  U  , 
since  l  >  2).  But  R/r  »  (R/p)(l  +  e  cos  v),  so  If  4s  -  a  -  a*  the  required  cancelling 
takca  place.  Thla  la  why  a'  is  the  preferred  1  . 

There  Is  no  corresponding  consideration  for  n  ,  and  the  obvious  course  is  to  identify  5 
with  the  secular  component  of  M  ;  thla  makes  it  only  weakly  dependent  on  the  definition  of 
M  itself,  as  M  is  independent  of  any  arbitrary  constant  in  M  .  We  work  from  (2),  but 


5 


nttd  to  eliminate  ch«  variability  of  n  ,  which  we  can  do  via  o’  ,  defined  auch  chat 
u'2a'^  -  u  •  Thua  (6>  gives 

n  ■  o'  -  3naD/y  ♦  O(U^)  ,  whence  ft  •  M  *  n*  +  9  -  3nSU/u  ,  (7) 

where,  aince  M  and  9  are  the  secular  componencs  of  Che  M  and  9  variation,  U  has  to 
be  interpreted  as  the  M-aean  value  of  U  .  To  derive  the  terms  of  (4).  we  develop  ft  , 

In  (7),  in  the  form  n’|l  +  At^l}  *  ^ut  before  pursuing  this  we  look  more 

closely  st  5  and  Z  for  the  ^-oaly  field. 

J2  analysis 

2 

Substituting  %<3f  sin  u  -  1)  for  P^Cain  6)  in  (6),  where  u  ■  v  +  w  ,  we  get 

2* 

»2  ‘  T5rf  /.  {if  <2h  +  3f  «*  2u>  «  •  (e> 

0 

Since  (p/r)^  dM/dv  -  ,  Integration  gives  U,  »  $y(K/p)$^h  to  first  order  in  J2  . 

Taking  9  -  Kftqh  /6/,  we  have  3nlU/u  ■  9  ,  so  (7)  gives  ft  •  n*  . 

Thus  the  preferred  values  for  ft  and  Z  are  n'  (for  t  *  2  only)  and  s'  (universally). 

2  3 

with  ft  Z  therefore  equal  to  unmodified  u  .  Kozal,  on  the  other  hand,  though  using  the 
same  mean  n  ,  has  a  different  mean  a  ,  denoted  here  by  ag  and  such  that 
4g  -  a’(l  -  -jR^h);  this  accounts  for  his  version  of  Kepler's  third  law.  As  Kozal  does  not 
give  s  rationale  for  eg  »  whilst  a  considerable  literature  follows  his  paper  / 4/  without 
clarifying  the  matter.  It  Is  worth  offering  an  explanation  here. 

The  attraction  of  an  unbiased  fir  was  noted  in  the  Introduction,  and  Kozal 'a  fir  is 
unbiased  for  circular  orbits.  Though  his  fir  expression  (for  general  e)  la  more 
complicated  than  the  one  to  be  given  here  (equation  (10)),  the  fir  of  (18)  is  biased,  even 

for  circular  orbits,  unless  h  -  0  (1  -  54.7*  or  125.3*).  The  lack  of  bias  If  ag  is  used 

means  chat,  for  an  orbit  that  Is  equatorial  as  well  es  circular,  r  (which  now  has  a  fixed 
value)  actually  Identifies  with  tg  (neglecting  J^).  The  value  of  a  (osculating,  and  also 

fixed)  Is  different,  however,  being  equal  to  Sq  ,  defined  In  the  next  paragraph.  Whilst 
considering  orbits  with  e  -  sin  1  •  0  ,  it  is  worth  remarking  on  a  third  candidate  (in 
addition  to  n'  and  tha  unbiased  n^)  for  ft  ,  viz  the  rate  of  change  of  longitude,  n^ 

say.  Since  din  (sign  given  by  cos  1)  Is  K5  ,  under  these  circumstances,  we  have 

*  o'  (1  ♦  K),  so  n^r  *11(1  +  K),  and  this  Is  of  particular  interest  for  a  geostationary 

orbit,  where  n.  Is  equsl  to  the  Earth's  rotation  rate;  if  a.  is  defined  such  that 

2  3  A 

n^  a^  -  u  ,  It  will  be  noted  that 

*0  -  •’  *  2<»'  -  V  *  2<‘x  •  V  •  (9> 

the  'unit  value'  of  jg  -  a^  being  about  0.52  ta  (geostationary  orbit  assumed). 

The  quantities  n^  and  Sq  ,  just  referred  to,  are  almply  the  M-mean  values  of  n  and  a  , 
but  were  introduced  (with  the  zero~sufflx  notation)  by  Kozal  (in  contradistinction  to  n' 
and  a^  ,  notated  as  ft  and  Z)  in  a  rather  misleading  way,  aince  he  describes  them  as 
'unperturbed  values  st  epoch'  though  they  are  actually  unbiased  values  that  are  Independent  of 
epoch  (the  word  'unperturbed*  seems  meaningless  In  the  context);  unfortunately,  this 
inapposite  description  is  followed  in  a  number  of  textbooks.  The  relations  between  (1) 
and  (Og  -  )  n',  and  (ii)  aQ>  Sg  and  s'  ,  are: 

(1)  oQ  -  n'(l-Kqh);  (li)  «0  -  ^(1  +  K5h)  -  .'(1+^)  .  (10) 

Analysis  for  (1  >  4) 

Cenersl  analysts  / 5 /  yields  expressions  for  9  (secular)  due  to  3^  ,  together  with  the 


6 


corresponding  ,  for  substitution  in  (7).  On  combining  the  two  terms,  for  each  t  ,  ve 

gee 

H  -  a'|1+^Jt(i)  Vlj  ’  vhere  vi  *  ;  (11) 

here  Aj(i)  is  *  polynomial  In  f  (if  1  is  even)  end  B^<e)  is  e  polynomial  In  e*  , 
normalized  such  that  A^(0)  -  B£(0)  *  1  ,  with  a  pure  number,  given  by  Cj  ■  \  end 

Cjj/C^  »  "  -(l  -  1)/1  if  t  >  4  .  (A^,  Bf  and  are  particular  cases  of  the  A*,  and 

j. 

Cj_  defined  in  /5/ ,  k  being  zero,  recurrence  relations  having  been  given  for  A  and  B  as 
well  as  C  .)  The  derivative  of  B^(e)  has  leading  term  fc(t  -  l)(t  -  2),  so  we  re-normallze 
(11),  writing  •  7153A1(i)Q1(i) ,  where 

Qt(«)  -  ^e_1dBl(e)/<ieJ /^(t  -  l)(i  -  2)j  and  If,'-  -  4(t  -  l)(t  -  2)C,  .  (12) 

The  recurrent.*  relation  for  Q£(e),  using  odd  values  of  1  as  well  as  even  ones,  is 

q,  -  [(21  -  5)0,^  -  (1  -  4)q2Qi_J  /  (t  -  l)  ,  (13) 

for  l  >  4  ,  with  Q2  *  0  and  Q3  -  1.  The  recurrence  relation  for  Y4  ,  using  even  values 
of  i  only,  is  (for  1  >  6  ,  with  •  9/8) 

7,  -  -[(»  -  l)2<t  -  2)/t(l  -  3)(l  -  *)  Yt.j]  .  (H) 

Since  Jjt  (for  t  >  4  and  even)  contributes  the  factor  1  +  2J£(R/p)*vj  to  (with 

S  -  a’),  we  can  write  (on  evaluating  the  y*  as  far  as  l  •  8) 


bV  ‘  v« 


.  735 
+  SZJi 


(j)9  Vd  -  "•  -  8?  *  s?2)j] 


-2 

(The  K  term  la  taken  from  /I / ,  where  it  is  not  derived.  A  full  derivation  should 
eventually  be  published  /&/•) 

3  MEAN  ELEMENTS  AND  COMPLETE  FIRST-ORDER  SHORT-PERIOD  PERTURBATIONS 

ASSOCIATED  WITH 

We  have  seen  that  the  expressions,  in  coordinates,  for  short-period  perturbations  depend  on 
the  definitions  of  mean  elements.  In  /6/  the  author  gave  expressions  for  the  first-order 
perturbations  due  to  an  arbitrary  harmonic  (tesserals  Included),  using  a  particular  system  of 
cylindrical-polar  coordinates;  complete  expressions  were  given  for  J2  perturbations,  but 
for  the  general  harmonic  (and,  additionally,  for  perturbations)  they  were  e-truncated, 

effectively  applying  only  to  near-circular  orbits.  More  recently,  the  results  of  a  complete 
analysis,  to  second  order  in  J2  and  first  order  in  J3  ,  have  been  summarized  /l/,  and  it  is 
Intended  that  full  details  will  follow  / 8/.  The  analysis  showed  that  a  spherical-polar 
system  is  actually  preferable  to  the  cylindrical-polar  system,  and  that  J2~based  definitions 
of  two  of  the  mean  elements  (3  and  8)  had  not  been  optimal  in  the  original  study  -  this 
conclusion  was  anticipated  by  the  'Note  added  in  proof  at  the  end  of  /6/.  Having  considered 
here  the  relation  of  the  optimal  1  and  ft  to  Kozal's  elements,  we  now  do  the  same  for  the 
remaining  elements,  concluding  with  the  resulting  expressions  for  first-order  perturbations 
In  coordinates. 

The  optimal  D,  0  and  M  are  identical  with  the  elements  of  Kozal,  so  we  only  have  to 
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consider  i  and  1  ,  but  It  is  slso  worth  covering  the  derived  element  $  .  In  terms  of  the 
optimsl  (preferred)  %  end  $  ,  Kozsi’s  elements  ere  given  by 

’  *[l  '  +  ?>]  “d  1>K  -  f[l  +  ~  3^)]  •  (16) 

The  rewerd  for.  using  our  preferred  1  (  ■  e*)  end  i  is  to  obteln  the  very  simple  expressions 
for  the  complete  perturbetlons  dr  end  6w  thet  follow.  Here  the  spherical-coordinate 
system  (r,  b,  w)  is  besed  on  the  instenteneous  position  of  the  mesn  orbltel  plane,  defined  by 

1  end  5  ,  such  thet  v  is  the  in-plane  eagle  (quasi-longitude)  end  b  is  the  out-of-plane 
engle  (quasi-latitude).  For  the  simplest  possible  expression  for  6b  (where  b  is  auto¬ 
matically  zero  by  definition  of  the  reference  plane) ,  the  1  we  need  is  within  0(K?)  of  the 
extreme  values  of  6  (at  the  apexes  of  the  orbit).  Kozal’s  inclination,  on  the  other  hand, 

_  2 

Is  within  0(KI  )  of  the  M-mean  1  .  This  seems  a  very  much  less  natural  choice  (and  is  a 
notable  example  of  the  allure  of  Che  intuitive  meaning  of  'mean');  it  differs  from  our 

preferred  1  by  an  0(K)  quantity,  since 

1^  -  I  +  kK  sin  2l  .  (17) 

The  mean  elements,  to  order  »  h*ve  now  »H  been  established,  and  the  complete  (first- 
order)  expressions  for  coordinate  perturbations  can  be  quoted.  They  are: 


6r  -  \  Kp(f  cos  2u  -  2h)  ,  <3b  •  f  Ke  sin  2i  { sln(u  +  ?)  -  3  sin  0  f 

and 

6w  »  TTK|f  sin  2u  +  4?f  sin(u  +  D)  +  8Sh  sin 


(18) 


The  corresponding  second-order  perturbations  are  not  given  here  (see  III)  but  it  is  worth 
noting  that  the  number  of  terms  required  for  6r,  6b  and  6v  are  6,  7  and  11  respectively. 
This  Is  in  striking  contrast  to  the  number  of  terms  required  for  the  six  osculating  elements, 
which  is  140  overall. 
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